TITLE OF INVENTION 

DERIVATION OF COMPOSITE STEP-FUNCTION RESPONSE 

BACKGROUND OF THE INVENTION 
5 The present invention relates to distance to fault measurements, and 

more particularly to the derivation of a composite step-function response from 
band-limited channel frequency response in distance to fault (DTF) 
measurements. 

For time domain reflectrometry (TDR) or distance to fault (DTF) 
10 transmission channel tests, step function testing along v^ith impulse testing is 
often used to measure wave propagation and reflections of the transmission 
channel. The step-function test is useful when the transmission discontinuity 
is frequency selective. Due to its time domain integration nature, the step 
function response is more sensitive to low frequency components. In 
1 5 measuring frequency response of a system the measured frequency range 
may be band limited, such as between 25 MHz and 2.5 GHz. 

Referring to Figs. 1a, 1b and 1c the respective graphs show the 
impulse response, the step function response, and a band limited step 
function response for a particular transmission channel. As readily seen, the 
20 step function response is very important in transmission channel diagnosis 

since the impulse response produces an almost insignificant difference where 
the reflection or discontinuity is a low frequency response, i.e., reflects low 
frequencies rather than high frequencies, while the step function produces a 
very noticeable difference as it contains mostly d.c. and low frequency 
25 components. 
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Generally the step function response may be derived by integrating the 
impulse response. As indicated above due to its integration nature, the step 
function response is more sensitive to low frequency components. For a 
frequency domain instrument, a transmission channel reflection coefficient is 

5 measured at each specified frequency within a specified range of 
frequencies, i.e., the measurement is band limited. TDR or DTF 
measurements are derived from the inverse Fourier transform (FFT'^) of the 
channel reflection coefficient response. When a frequency domain 
instrument cannot make measurements at low frequencies, i.e., the source 

10 covers a frequency range that excludes the low frequencies, an incorrect step 
function response may be produced when the band-limited TDR response is 
integrated, as shown in Fig. 1c, especially if the discontinuity or reflection is 
low frequency selective. 

What is desired is a method of deriving a composite step function 

15 response from a band-limited transmission channel response obtained from 
frequency domain measurements. 



BRIEF SUMMARY OF THE INVENTION 

Accordingly the present invention provides a method of deriving a 
20 composite step function response from a band-limited transmission channel 
frequency response. The method includes the steps of obtaining a time 
domain response from the band limited frequency response, identifying 
reflection events from the time domain response, estimating an impulse 
response from the identified reflection events, and determining the composite 
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step function from the estimated impulse response. The impulse response 
estimation is obtained from the observed time domain response as 

y(n) = h(n) - h(n)®w(n) 
where y(n) is the observed time domain response, h(n) is the impulse 
5 response to be estimated and w(n) is a window function 
w(n) = sin(a)o*n/Fs)/nn 
where coq is the initial frequency and Fg is the sample rate frequency. For 
reduction in calculation expense an impulse response segment is calculated 
over a narrow range of data about each reflection event. The resulting 
10 estimated impulse response is accumulated to produce the composite step 
response for the band limited transmission channel. 

The objects, advantages and other novel features of the present 
invention are apparent from the following detailed description when read in 
conjunction with the appended claims and attached drawing. 

15 

BRIEF DESCRIPTION OF THE SEVERAL VIEWS OF THE DRAWING 

Figs. 1a, 1b and 1c are graphical views respectively of (a) an impulse 

response, (b) a step function response and (c) a band limited step function 

response for a transmission channel according to the prior art. 
20 Fig. 2 is a graphical illustration view in the frequency domain of an 

estimated transmission channel impulse response algorithm according to the 

present invention. 

Fig. 3 is a graphical illustration view in the time domain of the 

estimated transmission channel impulse response algorithm according to the 
25 present invention. 
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Fig. 4 is a plan view of a display of the composite step function in the 
band limited transmission channel according to the present invention. 



DETAILED DESCRIPTION OF THE INVENTION 

5 Referring to Fig. 1c it is apparent that the band limited step function 

response does not give a true step function response to the discontinuity 
event or reflection that exists at that point. Therefore in order to provide an 
accurate step function it is necessary to estimate the transmission channel 
impulse response. The step function response derivation, as described 

10 below, has four steps: distance to fault derivation from band limited channel 
response; reflection surface detection or identification, impulse response 
estimation and step function response calculation. 

Distance to fault (DTF) is derived from reflection coefficients measured 
at discrete frequencies over a band limited frequency range. For 

15 computational efficiency and accuracy the following process is used, as an 
example: 

For B being the measurement bandwidth and F^ the upper frequency 
edge, then a working bandwidth B,, is 

U = IFJB] and B,, = F,/U 
20 A new central frequency becomes Fc = - 

and the lower frequency edge is F, = Fu - B^. 

Down shift the frequency band [F,, FJ to baseband (F,, = 0) and 
perform inverse discrete Fourier transform (IDFT or DFT"^) using sample rate 
F3 equal to 2B^. This is equivalent to directly performing IDFT on frequency 
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data. Then up-convert to frequency Fq = BJ2 and extract the real part of the 
signal as the DTF signal. 

In the digital domain this process may be simplified by letting H(a)k) be 
a reflection coefficient at frequency (k = 0, 1 , . . ., N-1), = 2nF| and co^.i 
5 = 2nFu. Then 

h{n) = (1/N)Ik=o.>N-iH(cOk)ei2"™" = IDFT(2H,2N) 
hpTF = Re(h(n)) 
H(a}o) = 0, . . ., H(a)^.i) = 0 where m={F^ - B)N/B^. 

The reflection surface identification may be automatic or user- 
10 interactive. In user-interactive mode a user inputs a center location (/| + /2)/2 
or the edges of the impulse response. In the automatic mode the center 
location is detected based on the event's reflection magnitude in the time 
domain. 

A detection function A(n) may be an envelope value of the up-converted 
15 signal /7(n) 

A(n) = \h{n)\ 

Alternatively using the local energy of the baseband signal as the detection 
function, 

x(n) = |/7(n)|sign(/7,(n)), A(n) = I„=.K.>Kx(n-m) 
20 where h{n) = h^{n) + jh^in). In either case if lA(n)| is greater than a threshold, 
then a reflection surface is detected and the center location determined. 

For impulse response estimation let h(n) denote the impulse response 
to be estimated and y(n) the corresponding band limited response (y(n) = 
hDTF(n))- Looking at Fig. 2 in the frequency domain Y(a)) = H(a)) - H(a))W(a)), 



where Y(a)) represents the band Nmited frequency coefficients, W((d) is a 
window function representing the frequency range not covered by the band 
limited source and H(a)) is the combination of the measured frequency 
coefficients and the frequency coefficients estimated over the window. In the 
time domain this becomes 

y(n) = h{n) - h(n)®w(n) = h(n) - I„h(m)w(n-m) 
the window function being w(n) = sinCcDon/FgVnn where Fs is the sampling 
frequency. As seen, the observed time response y(n) is a linear function of 
the impulse response h(n). Since y(n) may have many data points, it may be 
computationally expensive if every impulse response point is estimated. 
However normally each reflection surface covers a very short distance and its 
response lasts a limited time. Therefore for each identified reflection surface 
the impulse response is estimated over a narrow range of data around the 
reflection surface, as illustrated in Fig. 3. 

n^ < n < I, 

-Im=ii->iMh(m)w(n-m) 

is < n < nj 

y(n) = 

h(n) - Zm=ii.>iMh(m)w(n-m) i^ < n < I2 



In matrix form this is 

|y(ni) I |-w(ni-ii) -w(ni+1-ii) 



-w(-1) -w(-2) 
1-w(0) -w(-1) 
-w(1) 1d-w(0) 



1 
I 

I I 

|y(ii-i)l 
|y(ii) I 
|y(ii+i)l 
I ■ I 
i • I 
I • I 
|y(i2) I 
I • I 
1 • I 
I • I 

lyCns) 1 |-w(n2-ii) -w(n2-i,-1) 



-w(i2-ii) -w(i 



-w(ni+M-1-ii)| 
I 
I 
I 

-w{-M) I 
-w(1-M) ||h(mOI 
-w(2-M) \\h{m^)\ 
II ■ I 
II • I 
II • I 
1-w(0) ||h(mM)l 
I 



-w(n2-i2) I 



The iocalized h(m) may be optimally resolved by applying least-square error 
criteria. 

Y=DH, H = (D^D)-'D^Y 

The final step function response is calculated by accumulating the 
impulse response. 

x(n) = h(n)®u(n) = X„=_>^h(n-m)u(m) 
where u(n) is the step function 
1 n > 0 

u(n) = 

0 n < 0 

This produces 

x(n) = E„=o.>nh(n-m) 

Typical results of the band limited method described above are shown 
in Fig. 4 which shows the estimated impulse and step function responses of a 
band limited channel. The lighter line represents the estimated impulse 
response 20, showing an apparent small reflection 22 at approximately 200 
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meters and a larger negative impulse 24 at approximately 400 meters. The 
corresponding step function response 30 shows a significant response 32 at 
200 meters and a less significant response 34 at 400 meters. The 
combination of these responses gives a more complete picture of the events 

5 that occur in the transmission channel. 

Thus the present invention provides a method of deriving a composite 
step function response from a band limited transmission channel frequency 
response in distance to fault measurement by estimating the impulse 
response from the measured frequency response and a v^indow for the 

10 excluded frequencies, and accumulating the impulse response to obtain the 
step function response. 



